Oklo Constraint on the Time- Variability of the 
Fine-Structure Constant* 



Yasunori Fujii 

Advanced Research Institute for Science and Engineering, Waseda University, Tokyo, 169-8555 Japan 

Abstract 

The Oklo phenomenon, natural fission reactors which had taken place in Gabon about 2 bil- 
lion years ago, provides one of the most stringent constraints on the possible time- variability 
of the fine-structure constant a. We first review briefly what it is and how reliable it is in 
constraining a. We then compare the result with a more recent result on the nonzero change 
of a obtained from the observation of the QSO absorption lines. We suggest a possible way 
to make these results consistent with each other in terms of the behavior of a scalar field 
which is expected to be responsible for the acceleration of the universe. 



1 What is the Oklo phenomenon? 

Oklo is the name of the place of a uranium mine in Gabon, West Africa, near the equator. The 
mining company would supply the uranium ore to the French government. But in June of 1972, 
something unusual was noticed on the ore from Oklo; the abundance of 235 U was somewhat 
below the world standard, 0.7202%, well beyond the limit of permissible range. This might 
have undermined the company's reputation about the quality of their uranium. But finally 
after a few months of serious effort, French scientists came up with an unexpected, startling 
conclusion: 

The deficit of 235 U was a real effect of that a self-sustained fission reaction took place nat- 
urally in Oklo about 1.8 Gys ago, during the period of Proterozoic, part of Precambrian. In 
other words, natural reactors did exist well before 1942 when Enrico Fermi invented the artifi- 
cial reactor for the first time in Chicago. This has been called the "Oklo phenomenon," since 
then [1]. The result of their work was published in many ways, including [2]. 

See the separate figure, gabon2.jpg. 



Figure 1: Oklo in Gabon, West Africa, near the equator- 
There was a big press campaign, including the Le Monde article, for example, saying that 

Fermi was not an "innovator," but was only an "imitator" of Nature. 

Even more surprising was that the occurrence of this "natural reactors" had been predicted 

much earlier, 17 years earlier by a nuclear geochemist, Paul K. Kuroda in 1955 [3]. He discussed 
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several conditions. But the most important was that the ratio of 235 U, currently 0.70%, used to 
be much higher because of the different lifetimes of 235 U and 238 U; 7.038 x 10 8 y and 4.468 x 10 9 y, 
respectively. One can easily calculate the ratio 1.8 Gys ago to be as high as 3.2%. We note 
that 3% is a goal of most of today's enrichment facilities. Another condition was the presence 
of water which served as a moderator. 

So far 16 "reactor zones" (RZ) have been discovered in the Oklo area. In each of of them, 
extensive and detailed measurements have been made on the leftover fission products. 

2 How did Shlyakhter probe Aa? 

Under this circumstance, in 1976, Alex Shlyakhter [4,5] then in Leningrad proposed to look 
at 149 Sm, which is present naturally at the ratio 13.8%, but should be depleted in the reactor 
zones because it had absorbed neutrons strongly in the reactors 2 Gys ago, according to the 
reaction 

n+ 149 Sm^ 150 Sm + 7 . (1) 

One measures the abundance in Oklo reactor zones to estimate the cross section of this process, 
and compare the result with today's laboratory value. In this way one can tell how much 
nuclear physics 2 Gys ago could have been different from what it is. 

What is unique with this particular process (1) is that it is dominated by a resonance that 
lies as low as E r = 97.3mcV, while we know that a typical energy scale of nuclear physics is 
~ MeV. Compared with this, the above value is very small, nearly 7 orders of magnitude too 
small. This must be due to a nearly perfect cancellation between two effects; repulsive Coulomb 
force which is proportional to the fine-structure constant a, and attractive nuclear force which 
depends on the strong-interaction coupling constant squared a s . We are left with a very small 
leftover for the resonance energy, as illustrated Fig. 2. Suppose we change one of the coupling 
constants, a, say, only slightly. Then the strength of the Coulomb energy will change also 
slightly, and so will E r . However, the relative change may not be so small, because the starting 
value E r was already small. If this really happens, then the cross section may change rather 
significantly. This is a kind of amplification mechanism, which Shlyakhter exploited. 
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Figure 2: Nearly complete cancellation between the repulsive Coulomb energy and the attractive 
nuclear energy, leaving a small leftover for the resonance energy 

For the sake of illustration, we plot in Fig. 3 the cross section (based on the Breit-Wigner 
formula) as a function of AE r , the fictitious change of the resonance energy from today's value. 
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We also assume thermal equilibrium of the neutron flux; assuming a temperature corresponding 
to one of the curves shown. We find a sharp peak obviously coming from the resonance. Suppose 
E r at 2 Gys ago were smaller by lOmeV, a tiny amount. Suppose also T — 300°C, for example. 
Then we find the cross section bigger than today's value by a few %, a significant change. 

200 



* 150 



100 



50 



- n + u8 Sm 

: '/// \ 


1 II 1 1 1 II 1 1 1 II 1 1 II 

20°C 

100°C 

200°C 

300°0 

500°0 

1000"C 






- 


'•\\\ s • 



-200-150-100 -50 50 100 150 200 
Resonance Position Change AE r (meV) 



Figure 3: The cross section CT149 for the process (1) as a function of AE r — E r — 97.3meV, 
assuming thermal equilibrium of the neutron flux. The ranges of the observed cross section and 
the temperature are also shown, as given later by (6) and (7), respectively 



Incidentally, we are going to use the hat symbol attached to the cross section rather fre- 
quently. But this is only a technical convention of normalization, which is particularly popular 
among the researchers of the Oklo phenomenon. We do not worry too much at this moment. 

Taking advantage of this strong dependence of the cross section cr 14 g on AE r , Shlyakhter 
gave the upper bound 

\AE r \ < 50meV. (2) 

It is not very much clear how he derived this result, particularly how much the data uncertainties 
affected the conclusion. This is one of the points to be re-examined later. 

He still went on to discuss how this change of E r corresponds first to the change of the strong- 
interaction coupling constant, a s . He considered the resonance as a single-particle excitation 
in the potential, with its depth Vq ~ 50 MeV, which he assumed to be proportional to a s . If 
a s changes, Vq changes, and so does E r . Substituting from (2), he obtained the result 



Aa s 




AE r 


a s 




V 



50 meV 
50 MeV 



= 10" 



lcading to the value 10 9 . Further dividing by 2 x 10 9 y, he arrived at 

<5xl0~ 19 - 



y 



(3) 



(4) 



As for the electromagnetic interaction, he apparently replaced a s by a, resulting in dividing 
these by a/a s ~ 1/20, giving 

Aa 



a 



< 2 x 10" 



< 1 x 10~ 17 y~\ 



(5) 
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This last value has been several orders of magnitude more stringent than any other estimates, 
a kind of "champion result" for many years. However, one may raise a question against the 
argument from a s to a, and the suspicion may go further back to the derivation (3). 

For better understanding, we re-examined the whole analysis [6] by forming a team, which 
includes theorists of nuclear physics, reactor scientists and geologists. Among them, Hiroshi 
Hidaka is an expert nuclear chemist who has been specialized to the Oklo phenomenon. 

3 How good is it? 

We soon realized that a major error source of the data comes from the "post-reactor contam- 
ination," implying that certain amount of 149 Sm present in the outside environment flowed 
into the reactor core having occurred after the end of the reactor activity that is believed to 
have lasted several 10 5 years. This amount has nothing to do with what happened inside the 
reactors, so is a contamination from our purpose. This inflow was in fact the gradual mixture 
between inside and outside prompted by repeated successions of dissolution and precipitation 
of Sm, caused essentially by weathering, namely being exposed to the air. To minimize this 
embarrassing effect, we looked for samples in the reactor zones 10-16 discovered later than 
1984, deep underground, as shown in Fig. 4. Finally we decided to collect five samples taken 
from RZs 10 and 13 below the surface, with enough care of geologists expertise. 

See the separate figure, oklovrt2.jpg. 



Figure 4: Geological vertical cross section of the Oklo area 

The measured isotopic ratios related to 149 Sm for five samples are shown in Table 1. We see 
how small the abundances of 149 Sm are compared with the natural abundance 13.8%. We did 
not show the errors, but they are simply small. t\ is the time of the end of the reactor activity 
that started at t = 0. Also in the first line, we entered what is called "fluence," denoted by <f>ti, 
but actually a time-integrated neutron flux <j> during the whole duration of reactor activity. 

We then solved the evolution equations 

dN 147 (t)/dt = -ai47&N 147 (t) + N$ 35 exp(-a a (f>t)af235(i>Yi47, 
dN U8 (t)/dt = u U7 ^N U7 {t), 

dNu 9 (t)/dt = -<Ti4907Vi4 9 (t) + A^ 35 exp(-CT a 0t)o-f2350Y"l49- 

for the related isotopes to calculate the cross section aug for the process (1). The result is 
summarized, 

(7149 = (91 ± 6) kb, (6) 
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Table 1: Measured isotopic ratios related to 9 Sm obtained from five samples. The fluence and 
cross section arc also listed 



Sample 





SF84-1469 


SF84-1480 


SF84-1485 


SF84-1492 


SD.37 


hx (Vkb) 


0.525 


0.798 


0.622 


0.564 


0.780 


^144 (*l) (%) 


0.1052 


0.2401 


0.2073 


0.1619 


0.06909 


N U7 (h) (%) 


55.34 


53.23 


54.03 


54.81 


52.74 


Nus(h) (%) 


2.796 


3.468 


3.079 


2.890 


4.694 


^149 (ti) (%) 


0.5544 


0.2821 


0.4466 


0.4296 


0.3088 


JV236(*l)/JV238(tl) 


0.03181 


0.02665 


0.02971 


0.03047 


0.02435 


(Tug (kb) 


85.6 


96.5 


83.8 


99.0 


89.5 



corresponding to the narrow horizontal band shown in Fig. 3. 

We also made an estimate of the temperature, by the traditional way supplemented by a 
latest technique, giving 

T = (200 - 400)°C, (7) 

corresponding to the shaded area in Fig. 3. We find two intersections, and the corresponding 
two separated ranges of AE r . 

( (9 ± 11) meV, right-branch, Null 
AE r = { (8) 
{ (-97 ± 8) mcV, left-branch, Non-Null 

The right-branch range covers zero, so that a null result in the usual sense, while the other 
implies that E r was different from today's value by more than 10 standard deviations. Does 
this really imply an evidence of the difference in 2 billion years ago? We tried to see if the 
non-null result can be eliminated by looking at other isotopes like 155 ' 157 Gd, but so far no final 
conclusion yet. 

At this point we compare our result with those due to Damour and Dyson [7] (DD), who 
used the samples obtained mainly from near the surface, giving the cross section: 

ai4 9 = (75 ± 9)kb, (9) 

somewhat smaller than our result (6). This seems consistent with our suspicion that their data 
suffered from contamination. Also, they did not come to separating the two ranges. They could 
have done it, though the "right-branch" range would failed to cover zero even at the level of 2 
standard deviations. Instead, they gave only a combined range, —120 meV < AE r < 90 meV, 
which more than covers our two ranges. 

We admit that we are still short of determining which range is correct. Then one might say 
that we should also be satisfied by the combined range. But we still insist that it is a progress 
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to have established a disallowed range in between. 

Now we move on to discuss how AE r is translated into Aa following DD. First they ignored 
the contribution from the strong interaction entirely, focusing on the first "term" coming from 
the Coulomb contribution in Fig. 2. 

Consider the energy A4 C , given by the difference of the Coulomb energies between the states 
with 150 and 149. They paid special attention to the fact that the resonance in the 150 Sm is 
excited. But we simplify the analysis, at this moment, by appealing to the semi-empirical mass 
formula due to Weizsacker, finding M. c « — l.lMeV. In this calculation, one has to allow an 
error perhaps within the factor 2 or 1/2. Notice also that the above result is negative, apparently 
in contradiction to the illustration in Fig. 2. Obviously we dealt basically with a repulsive force, 
but we calculated the difference, which turns out to be negative. Nothing is wrong, but we would 
better put the two terms upside-down on the left-hand side, but keeping the right-hand side 
still positive. We may also assume that M c oc a. We then obtain AE r = AM C = (Aa/a)M c , 
thus giving 

Aa _ AE r _ ( (-0.8 ± 1.0) x 10~ 8 , null, upper bound 
a M c ~ \ (0.88 ± 0.07) x 10" 7 , non-null 

Divide by -2 x 10 9 y to get 

a ( (0.4 ± 0.5) x 10 _17 y _1 , null, upper bound 
a = \ (-0.44 ± 0.04) x 10~ 16 y~\ non-null 

This upper bound happens to agree quite well with Shlyakhter's result 1 x 10 _17 y _1 . The 
agreement to this extent seems, however, rather accidental, because, among other things, it is 
unlikely that the data as good as ours was available in 1976. 

We emphasize here that the simple estimate due to DD, as described here, might be called 
"Coulomb-only estimate," which serves as a basis for more general analyses. 

In fact what really happens might be a combined result of both interactions, and one wishes 
if one could include the strong interaction as well. But then everything is going to be compli- 
cated, for example, like the QCD analysis in [8] . But there are something independent of such 
complications as long as we appreciate the condition that AE r is much smaller than either of 
the mass scales, M c and its strong-interaction counterpart M s - First we find from Fig. 2 that 
the mass scale of M. s is nearly equal in its size to that of A4 C , obviously much smaller than 
Shlyakhter's "50 MeV." It then also follows that Aa s /a s should be nearly of the same size as 
Aa/a. Of course there arc some differences from the Coulomb case; M s may not be simply 
proportional to a s . This may result in a revision of a factor, but certainly not of an order of 
magnitude. Then we go through a bit of analysis to conclude finally that it is unlikely that, by 
the strong interaction, Aa/a deviates from the Coulomb-only estimate by more than an order 
of magnitude, no matter how complicated the exact analyses might be. It can be smaller. See 
Appendix A for more details. 

Then, as always, there is a possibility of an exception, no matter how remote. This allows, 
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in principle, that both of Aa/a and Aa s /a s are quite large, in fact without limit, but cancel 
each other leaving a small value of AE r . At this moment, however, we assume that no such 
fine-tuning nor coincidence occurs in the real world. 

Let us summarize what the situation is with the question of "uncertainties." In the theoret- 
ical aspect as a whole, we say again that the Coulomb-only estimate of the relative change of 
a is correct likely within an order of magnitude, for whatever the complication of the effect of 
strong interaction. We only add a few related remarks. 

• Only AE r and the related cross section aug are sensitive to Aa. No other quantities are. 

• Our formulation is such that, to a good approximation in practice, the neutron flux can 
be any function of time. For example, the reactor activity can be even "sporadic." The 
only thing that counts is the fluence, no matter how long it took. 

• Estimating fluence is complicated but is a standard estimate and is reliable. 

• In principle we may not rule out one of the higher resonances to come down near the 
threshold, giving much larger value of AE r . We have a reason, however, to believe this 
to be highly unlikely. See Appendix B. 

In the observational aspect, we repeat our previous statement; post-reactor contamination 
is the largest error source, with a few more comments; 

• Nothing is serious for 149 Sm in our samples. (A few % contamination seems even better 
because then the range of AE r for the null result covers zero more in the middle.) 

• This is not the case for 155 ' 157 Gd which enjoy even lower resonance energies. The ab- 
sorption cross sections are also larger. But the effects are too strong to the extent that 
the residual abundances are too little, so are too sensitive to contamination, even with 
our samples with minimized effect of weathering. Shlyakhter was clever, when he chose 
149 Sm. This is also precisely why we reached short of complete elimination of the non-null 
result, as mentioned before. 

4 How can it be consistent with the QSO result? 

According to V. Flambaum and M. Murphy at this meeting, the latest version of their result 
on the time-variation of a from spectroscopy of QSO absorption systems is [9] (see also their 
contribution to the proceedings [10]): 

Aa 

= (-0.54 ±0.12) x 10~ 5 . 

a 

We show in Fig. 5, taken from Fig. 8 of [9], our own plot as a function of the fractional look- 
back time u defined by u = (t — t)/t , with t the present age of the universe. Their weighted 
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mean can be viewed as a fit by a horizontal straight line at —0.54, as also shown in Fig. 5. We 
will call this a "1-parameter fit" for the later convenience. Notice that x 2 cd = 1.06. 




Figure 5: QSO result from the 128 data points is shown in the upper panel, while the 13 
binned data might provide an easier view in the lower panel [9] . The long-dashed line is for the 
weighted mean —0.54. The Oklo time u kio = 0.142 is also shown 

We now include the data from the Oklo phenomenon, at u kio = 0.142, barely outside the 
QSO range, as also shown in Fig. 5. In this sense they are different things. However, we should 
put error-bars, which are invisibly small in this plot; less than 10 -4 and 10~ 2 — 10~ 3 in the 
horizontal and vertical directions, respectively. Remember that Aa/a from the Oklo constraint 
is 10~ 7 — 10~ 8 . If we extend the straight line naively down to the Oklo time, it will miss the 
point off (10 2 — 10 3 ) standard deviations, resulting in an enormous value of x 2 , too large to 
be acceptable. One wants to bend the line to pass the point almost exactly, but one needs a 
physical reason. What is that? This is the issue. Already there have been several attempts 
[11-16]. But we stick to our own idea that this issue has something to do with the accelerating 
universe, another big issue in today's cosmology. 

Now probably everyone knows that our universe is accelerating [17]. This behavior is best de- 
scribed in terms of a positive cosmological constant, whose size is given usually by the parameter 
Oa = A//9 cr ~ 0.7, where the critical density p cr is given by ~ t^ 2 . The coefficient here is of the 
order one if we use the reduced Planckian unit system in which c = h = Mp(= ^/87rG/ (ch)) = 1. 
In this unit, the present age to ~ 1.4 x 10 10 y is about 10 60 . So we find p cr ~ A ~ 10~ 120 . 

Today's cosmological constant problem has two faces or questions: Why is it so small? Why 
is it still nonzero? The first question can be replied by the "scenario of a decaying cosmological 
constant;" A is not a true constant but decays like ~ t~ 2 [18,19]. This simple idea can be 
implemented by the "scalar-tensor theory." We expect that a scalar field plays an important 
role. This scalar field may have its origin in string theory in which a graviton has a spinless 
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companion called the "dilaton." 

The second question seems to require a deviation from the simplest version of the scalar- 
tensor theory. As one of the possible ways we call for another scalar field, called \, in addition 
to the dilaton a. These two fields comprise what is called the "dark energy," and their energy 
density p s — p(a,x) is interpreted as an effective cosmological constant A ff. There are many 
details involved, though we are not going into any details. Readers are advised to refer to our 
recent book [20]. 




20 40 

log t 



Figure 6: An example of our cosmological solutions in the Friedmann universe wit k = 0. The 
scale factor a, the scalar fields er and \i the energy densities p and p s of the ordinary matter and 
the scalar fields, respectively, are plotted against \ogt, in the reduced Planckian unit system. 
The present time is around logi w 60. The middle panel shows an effective exponent of the 
scale factor, In a/ In t 

Skipping all the details, we show in Fig. 6 an example of our solutions in the Friedmann 
universe with flat 3-space. The horizontal axis is logt. In the Planckian unit system, the present 
time is somewhere around 60. In the lower panel, p is the usual matter energy density, which 
falls off roughly as t~ 2 . The energy density of the scalar fields p s is the effective cosmological 
constant, A e ff, also falling off like t~ 2 as an overall behavior, thus respecting the scenario of 
a decaying cosmological constant. 1. But the plot also shows occasional deviations, notably the 
plateau behaviors. Obviously each plateau mimics a cosmological "constant." Furthermore, it 
comes to a crossing with the ordinary matter energy density. One of them is expected to occur 
around the present epoch. Nearly in coincidence with this crossing, we find a "mini-inflation" 
of the scale factor; a bit of sharp increase in In a shown in the upper panel. This nicely fits the 
observed acceleration of the universe. 

tAs emphasized in [20], the gravitational "constant" G in the physical conformal frame identified (nearly) 
with the Einstein frame in the this model is (nearly) time-independent, instead of decaying with time as in [19] 
presented in the Jordan frame 
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Also shown in the upper panel are the sudden changes of the scalar fields a and x, again 
in coincidence with the crossing between p and p s . However, the most interesting is to take 
a close-up view of what appears to be a simple and small jump of a. With the magnification 
rate as large as 330 in the vertical direction, we find a surprising behavior shown in Fig. 7, 
something like a damped oscillation. 




Figure 7: A close-up view of a and \ m the upper panel of Fig. 6 around the present epoch. 
The magnification rate in the vertical direction is 330 

We know what the underlying mechanism is at the deeper level. But we still say that this 
damped-oscillation-like behavior is in fact "the heart and soul" of the entire dynamics that 
eventually brings about the acceleration as we see it. On the other hand, the acceleration itself 
docs not care how invisibly tiny oscillation is taking place behind the scene. There are many 
variations in the way of oscillation. In this sense, we have a degeneracy, which the cosmological 
acceleration does not resolve. However, this invisibly small oscillation may show up through 
the time- variation of a. 

This is an expectation based on a general view that changing a if any is due to the changing 
scalar field, expressed symbolically as 

— cx Act. (10) 

a 

String theory suggests this dependence for the gauge coupling constant. We ourselves derived 
a relation of this type, based on QED, featuring a quantum-anomaly type of calculation. But 
we do not want to be too specific on these theoretical details, nor to depend heavily on the 
choice of the solutions, like the one in Fig. 7. This is particularly crucial because we have many 
different solutions for a given cosmological behavior we want to fit. Rather, we are going to 
follow a phcnomenological approach which we describe briefly, leaving more details to [21]. 
Let y denote Aa/a in units of 10~ 5 . Then we assume a dependence on the fractional 
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look-back time u in the way of a damped oscillation 

y = ae ou sin I 2n J, (11) 

where the parameters are going to be determined to fit the QSO data as well as the Oklo 
constraint. 

We first choose M kio = 0.142 corresponding to the Oklo time of 1.95 Gys ago. The Oklo 
constraint, to be 10~ 2 — 1CP 3 in terms of y, is approximately zero in this scale. The remaining 
parameters a,b,T are determined by minimizing x 2 f° r the QSO data. In this sense we call this 
a 3-paramctcr fit. We do not include the Oklo data in computing % 2 , because we consider the 
Oklo has been already fitted approximately by choosing a zero of the function as above. 




0.2 0.4 0.6 0.8 

Fractional look-back time u 



Figure 8: The solid curves are for the 3-parameter fit with the least x 2 ( a = 0.151, b = 2.4, T — 
0.714, x 2 e d = 1-09), to be compared with the 1-parameter fit, represented by a horizontal 
straight (long-dashed) line at —0.54 in Fig. 5. Note the 10 times magnification for u < 0.2. 
The dotted and dashed curves are for b = 0.0 and b = 4.0, respectively 

We limit ourselves to a region of a, b, T in a manner roughly consistent with the theoretical 
model of the accelerating universe. In this range we searched for local minima of x 2 . Among 
several of them we find the least minimum which is given by b = 2.4, a = 0.151, T = 0.714 
resulting in x 2 c d = 1-09- This x 2 c d i s similar to x 2 c d = 1-06 obtained for their 1-parameter 
fit. In this sense our 3-parameter fit is nearly as good as the fit in [9]. The solid curve in 
Fig. 8 shows the actual plot, probably better shown in the binned plot in the lower panel. We 
magnified the curves below u = 0.2 by 10 times. 

One might ask us why we are satisfied with x 2 cd = 1.09 which is not smaller than 1.06 for 
the 1-parameter fit, in spite of the fact that we have more degrees of freedom. We answer the 
question by pointing out the following: 

• Our Xred = 1-09 is for the whole data including the Oklo, because, as we noted, the 
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Oklo constraint has been already "included" in a sense. In comparison, however, the 
1-parameter fit gives an unacceptably large x 2 when we include the Oklo, which was the 
starting point of the whole discussion. 

• Our 3-parameter fit was motivated originally by a theoretical "prejudice." There was no 
guarantee that it fits the reality We are relieved to find that our prejudice somehow 
survived a realistic test. 
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Figure 9: The 3-dimensional 68% confidence region is illustrated in terms of three cross sections 
for b = 0.0, 2.4, 4.0, as marked beside each contour, shown in 2-dimensional a — T space 

We now discuss what the confidence region for 68% is like for this 3-parameter fit. We 
imagine a confidence volume in 3-dimensional space of 6, a, T. Figure 9 shows, however, 2- 
dimcnsional cross sections for three different values of b. The contour labeled by 2.4 in this 
figure shows the one for b — 2.4. We also show other 2 cross sections corresponding to b = 0.0 
and 6 = 4.0, respectively. They give only slightly larger \ 2 than that for b = 2.4. In this way 
one imagines what the 3-dimensional volume looks like. The curves for these b are also plotted 
by the dotted and dashed curves, respectively, in Fig. 8. They are different from each other 
only in the lower- u region, u <, 0.5. 

We further add that we obtained several other solutions with other values of b, a, T which 
give local minima of x 2 , as we indicated before. As it turned out, however, they tend to give 
This is a number which is nearly comparable with x 2 od = 1.24, which we would 
obtain by fitting the QSO data by a horizontal straight line y = 0, namely the w-axis itself. We 
may have a good reason to exclude these fits. 

We may compare the result shown in Fig. 8 with the cosmological solution which we started 
from, as we showed in Fig. 7. For the latter we may estimate the parameters approximately, 
which will be shown in Table 2, together with the corresponding ones for the former. 

Agreement in the values of b is obvious. We have to have more theoretical details in order 
for the comparison of a to make sense, though we may reasonably find a consistency. On the 
other hand, there is a discrepancy between the values of T. The available QSO data shows a 
rather flat distribution of Aa/a, which favors a "larger" T. In this connection we point out, 
however, that we have chosen the solution in Fig. 7 rather arbitrarily. In fact Fig. 10 indicates 



b 


a 


T 


X s 


A. 


0.0 


.728 


.716 


136.83 


1.10 


2.4 


.151 


.714 


136.19 


1.09 


4.0 


.051 


.71 


136.46 


1.09 



"4.0 2.4 o.O 
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Table 2: Comparison between the damped-oscillator-like fit with an example of the cosmological 
solution 



a 



Cosmological solution ~ 2 4 ~ 2 5 ~ 22 
in Fig. 7 

Fit in Fig. 8 0.15 ±0.05 2.4 ± 2 0.71 ± 0.06 



that it happened that we have come across a relatively small T. We only conclude that we 
have to look for other solutions of the cosmological equations which still fit the way of the 
cosmological acceleration. 
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Figure 10: Another magnified view of the behavior of a. The previous Fig. 7 is only a small 
part of the present figure, corresponding to a curve with the attached number 7544. These 
numbers show the last 4 digits in the initial value at t = 10 10 , as explained in detail in [20]. 
The nucleosynthesis and CMB times are « 45 and « 55, respectively. The vertical scale may 
be accepted arbitrary at this moment 

Finally, we add the following comments for further improvements of the fit. 

• We show Fig. 10, taken from Chapter 5 of [20], which includes Fig. 7 as a small portion 
for a special choice of the initial value. This figure also demonstrates that the detailed 
behavior of the small oscillation depends heavily on the early history, depending sensitively 
on the behavior of the scalar fields at the initial times, particularly on those around the 
time of nucleosynthesis. In this sense determining Aa/a in this epoch as well as in the 
CMB era is crucially important in this approach. 

• In our approach in terms of (11) we assumed that the present time corresponds to the 
limiting, still transient, behavior toward the common flat value immediately prior to the 
big and sudden jump of a. It seems better to consider that we are in the middle of the 
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oscillation behavior, in general. Taking this possibility into account will make it more 
likely to satisfy the natural condition y(0) = 0, which should be true by definition. This 
type of the fit will be discussed in Appendix C. 

• Our analysis is based on the simplest assumption on the a — \ interaction, as given by 
(5.58) of [20]. This might be modified to improve the fit. 

As the very last comment, we hope if natural reactors will be discovered somewhere else, 
thus providing us with additional constraints, hopefully at different times. 

Appendix 

A Bound on Aq/q from the Coulomb-only estimate 

The situation described in Fig. 2 may be given the expression 

E r =M C + Ms, (A.l) 



with the condition 



We then obtain 



a a s 
In deriving the second equation we assumed 



\E r \ « \M c (a)\ ~ \M s (a s )\. (A.2) 



9E r . dE r . 
AE r = —Aa + —Aa s , 
oa oa s 

^ C -Aa + — Aa s . (A.3) 



dM c M c 8M S M s 

— — = , and — = , (A.4) 

oa a oa s a s 



to simplify the equations, for the moment. 

According to (A.l) and (A.2) we put (A.3) into 



a a 

Ignoring the second term yields the Coulomb-only estimate 

Aa AE, 



(A.5) 



D c0 . (A.6) 



a M c 

On the other hand, we notice that the right-hand side of (A.5) happens to vanish if 

Aa s _ Aa 
a, a ' 



(A.7) 
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leaving Aa/a undetermined in terms of AE r . (We then have to bring E r back again on the 
right-hand side of (A. 5), as discussed at the end of Section 5 of [6], thus corresponding to the 
"exception" mentioned toward the end of Section 3. We ignore this case at this moment.) We 
may assume, however, 

^=^, (A.8) 

a s a 

where £ is to be determined based on the more fundamental laws of physics, as attempted in 
[8,22,23]. By using this in (A. 5), we obtain 

— «(l-0 -1 A*. (A.9) 
a 

Suppose a special relation (A. 7) holds true within the accuracy of 10%, for example. This 
implies that (A.8) holds true for 

£=l + <5£, with < 0.1. (A.10) 

Then (A.9) implies 



1 Aa 

D c o a 



|^| _1 >10. (A.ll) 



In other words, Aa/a should remain close to the Coulomb-only estimate within an order of 
magnitude, unless the equality Aa/a = Aa s /a s holds true to the accuracy better than 10%. 

This result may be extended to more general situations, in which A4 C and A4 S depend on 
a s and a, respectively, though then separating into A4 C and A4 S in (A.l) may not be unique. A 
certain relation like (A. 7) is expected to result in the vanishing right-hand side of an equation 
corresponding to (A. 5). It is unlikely that the relation of this kind holds true exactly in practice. * 
Unless it does within the accuracy of 10%, we should always expect Aa/a to remain less than 
an order of magnitude of the Coulomb-only estimate. 

B Distant migration of the higher resonances 

We have so far assumed that AE r is very small, much smaller than E r = 93.7 meV for Sm. In 
fact we obtained AE r as small as 10 meV, thus giving Aa/a ~ 10~ 8 . This is a right attitude 
as long as we try to find as small an upper bound as we can. Now, however, the QSO result 
indicates a much larger value, up to 2 or 3 orders of magnitude larger. This might raise a 
question if the Oklo phenomenon does in fact yield a correspondingly larger AE r . One may 
suggest that we are looking at in the remnants of Oklo RZ a distant "migration" of a higher 
resonance down nearly to the threshold of n + 149 Sm. 

This possibility was already discussed by Shlyakhter [5], based on a statistical argument. 
Inspired by his approach, Akira Iwamoto and the present author attempt here a similar analysis 

■fProm the minimal supersymmetric standard model follows £ ps 6 at fi = Mz [23], though this result is not 
readily extrapolated to a much smaller fi in QCD, making it even unlikely to derive a value anywhere near unity 
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on the issue, reaching a rather negative conclusion by including the observation of Gd. We 
start, however, with discussing Sm first, for which the energies and the widths of the first four 
resonances are shown in Table 3. 

Table 3: First four resonances in n + 149 Sm. T n is the elastic width. The last line represents 
approximately expected time-variability of a obtained by the Coulomb-only estimate in units 
of 1(T 5 

Resonances 1st 2nd 3rd 4th 

E r (mcV) 97 872 4950 6430 

T„ (meV) 0.53 0.74 2.13 0.72 
r tot (meV) 30 30 31 30 
Aa/a (10~ 5 ) - 10~ 3 0.08 0.45 0.58 



Suppose AE r was negatively so large 2 Gys ago that one of the higher levels came down 
nearly at the same position as the first resonance. We may expect a significant effect if its 
(shifted) energy minus 97 meV happens to be close to the threshold. We may ignore the 
thermal energy (60 — 90) meV corresponding to the estimated temperature of (200 — 400)°C in 
the natural reactors. According to Table 3, the elastic width T n is found to be proportional 
to the center-of-mass momentum within an order of magnitude. We then expect the Coulomb- 
only estimate Aa/a = AE r /M. c ~ —E r /M c to be roughly correct with a common value 
M c « —1.1 McV for E r <C McV. This is the way we have obtained the last line of Table 3. 
Interesting enough, the values for the last two resonances turn out to be comparable with those 
reported by the QSO observations [9], but with the wrong sign. 

We notice, however, that r to t is considerably smaller than the energy difference required 
for the shift, implying that the "probability" of finding a shifted level that falls in the range of 
~ Ttot ~ 30meV around the threshold is rather small. This conclusion seems to be corroborated 
by taking 155 ' 157 Gd into account as will be discussed. 

There are many excited levels also in n + 155,157 Gd, some of which are illustrated in Fig. 
11. The first resonance levels appear at exceptionally low energies, 26.8meV and 31.4mcV, 
respectively. They are nearly degenerate. This unique feature is shared by none of the higher 
resonances, though all of them show remarkably similar widths, around 100 me V. 

The analysis of the Oklo natural reactors shows, however, that considerable enhancement 
near the threshold appears to occur for both isotopes. This conclusion seems to remain true even 
if possible significant effect of "contamination" is included [6]. What we observe is reasonably 
interpreted as coming from both of the two resonances. 

The inherent ambiguity coming from higher resonances as encountered for n + 149 Sm is 
also relevant here. This time, however, we must expect the levels to land near the thresholds 
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Figure 11: Levels of resonances of n+ 15 ^3d and n+ 157 Gd, shown above and below the horizontal 
axis, respectively. The scale of the energy levels is eV. They appear to be distributed rather 
randomly except for the first resonances which are nearly degenerate 

"simultaneously" in both reactions. This requires "squared" smallness of the probability in 
view of Fig. 11. Note that the ratios of the widths to the level spacings are even smaller than 
those in Sm. Small probability will be "cubic" if we combine all of the results of Sm and Gd. 

All in all, it is highly unlikely that the "discrepancy" between the Oklo constraint and the 
QSO result can be removed by assuming "distant migration" of higher resonance levels of the 
relevant isotopes. 

C Another 3-parameter fit with an offset 

We have recently found a fit with an offset parameter included, parametrized by 

y(u) — a (e bu cos (v — vi) — cos («i)) , (CI) 

where v/u = w kio/w kio = with v\ determined by 

wi = tan" 1 (( e - fe «okio _ cos(w ok io)) / sin(w oklo )) .Is (C.2) 

We easily find that y defined this way vanishes both at u = and u = w kio- The three 
parameters a, 6, T are then determined to minimize \ 2 f° r the QSO data. The result is for 
a = 0.046, b = 4.0, T = 1.307 with the fit shown in Fig. 12. The resulting x r 2 cd = 1-071 is even 
smaller than 1.09 for our previous 3-parameter fit in Fig. 8. 

Through these analyses we also find it unlikely that the current QSO result respects another 
constraint at the "meteorite time" around 4.6 Gys ago, or u « 0.33, requiring \y\ < 0.025 [8]. 
We are re-examining the basic formulation used in analyzing the decay of 187 Re ^ 187 0s [24,25]. 
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Figure 12: The best fit of (C.l) for a = 0.046,6 = 4.0, T = 1.307 giving x? od = 1-071 which is 
comparable with 1.06 for the weighted mean of [9] 
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